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AB.9TRACT  -  Often  it  t«  possible  to  decompose  a  wave  equation 
into  vertical  and  borizoatal  component*.  It  b  osefal  to  cooler 
Ibe  case  when  the  vertical  dependence  b  completely  expressed  ta 
a  vertical  normal  mode  representation.  A  parabr^ie  wave  cquatiMi 
with  ran^e  dependence  transforms  to  6rst  order  linear  diiTereolUl 
eqoations  in  an  inBnite  dimensional  Hilbert  space.  Tfab  permits 
concrete  expression  of  complicaled  functions  of  operators  and  thdr 
calculation.  Tbe  non-constant  coefScient  matrix  shows  range  de¬ 
pendent  vertical  mode  interactions  and  constraints  m  amplitudes 
In  thb  study  difTerences  and  sirrulariUes  between  sarious  trans¬ 
formed  equations  ate  explored,  lie  proposed  method  also  permits 
comparison  to  exbting  Hilbert  space  aoal>sb  of  the  system  when 
the  sjslem  b  slq^astk. 

1.  INTRODUCTION 

Many  authors  formally  derive  various  parabolic  equations  in 
the  context  of  approximaUng  the  Helmholtz  equation  in  ocean 
acoustics  (see,  e  y.,  McDaniel  (I]  or  Collins  (2]).  One  can,  however, 
use  some  function  analytic  methods  in  Sturm-Liouville  theory  to 
make  explicit  theoretical  calculations  (Keller  and  Ahluwalia  (3]). 
McDaniel  has  used  some  normal  mode  theory  to  discuss  various  er¬ 
rors  of  the  parabolic  approximation  [6].  Furthermore,  calculations 
with  (he  generahred  PourKr  transform  gtte  concrete  realtzations  of 
functions  of  operators.  These  allow  the  statement  of  cetlam  tech¬ 
nical  conditions  on  the  range  of  applicability.  Although  we  start 
with  Tapperl's  equation  (f],  the  method  works  for  many  split  step 
methods  and,  with  certain  timilalions,  (he  eilipUe  wave  equation. 
Wt  also  restrict  our  attention  to  a  flat  (op  and  bottom  surfaces 
with  constant  Sturm-Ltouville  boundary  conditions.  These  approx¬ 
imate  pressure  release,  rigid  bottom,  or  (he  Peketis  (wx)  fluid  la)er 
trapped  modes  (e/.  Kinsler,  Frey,  Coppens,  and  Sanders  (3))  Thb 
will  lead  to  a  completely  discrete  spectrum,  i  e ,  trapped  modes 
cnt),  but  the  procedure  b  exactly  analogous  when  there  b  contin¬ 
uous  spectrum  corresponding  to  dtlTracled  modes  The  latter  lead 
to  integral  operators  rather  than  matrices,  but  with  the  appropri¬ 
ate  changes.  These  differences  make  for  slightly  dilTerent  technical 
conditions  from  (he  discrete  case. 

?.  THANSTORMING  TUB  PARAIIOf.lC  EQUATION 

We  start  with  Tappetl's  original  parabolic  equation  (I) 


Since  the  complete,  we  may  write  the  generalized 

Fourier  series  foe  u  and 

The  inner  product  representation  for  the  Fourier  cocfHcienl  is 

bi„  =  hf^{r)  =  {u,^„)  =  y  ^«(r,z)^„dr.  (6) 

IVansformatbn  of  the  parabolic  equation  requires  r  derivatiies,  so 
differentiate  (3),  to  And 

+  +  +  =0  P) 

Now  use  the  series  for  (5b)  m  Ih^  square  brackets 

+  («  +  +  (fl  r  +  ^j=0 

Rewrite  the  dj  -f  f  m  the  bracket*  as  an  eigi^malue  using  (3)  F< 
nall>,  (he  orlhonormalil>  of  the  eigenfunctions,  ,dn)  ~ 
allows  us  to  pull  off  individual  terms  in  the  senes.  This  results  m 


(  0,  jszp 

(6) 

(9), 

We  calculate  0  s  =  2o^.  As  a  special 

case,  suppose  the  (he  change  m  potential  y  r  »  independent  of  z, 
t  e.  s  0,  then  aj^  s  0  for  all  y  and  p,  and  if,,,  s  for  all  p 
We  now  combine  the  previous  compulations  to  And  the  generalized  » 
Fourier  transformed  difTerential  operations  on  u:  I 

in  (he  r  direction,  and  in  the  z  direction 


(2iI5.  +  5j  +  ?)u=:0  (I) 

uheref  =  =  I*(n*  —  1)  has  r  dependence  (ignoring ap¬ 

proximation  to  (he  elliptic  Ilelmhollz  equation  initially).  We  write 
d,/  s  The  boundary  conditions  will  be  r  independent 

cosf?(  u(r,<) -sin9^  u,,(r,<)  s  0  (2) 

where  C  =  O.^h  =  0  OR  a  pressure  release  surface).  For  a  ftied 
r,  one  may  consider  (he  Stuinvhiouville  problem  in  the  z  variable 
(see,  c  s  >  Coddinglon  and  Levinson  (7]  or  Weidmann  (8)) 

+  +  »/<>/  =  0  (3) 

where  salisAes  the  same  boundary  conditions  (2).  For  continu¬ 
ous,  relatively  smalt  g  one  may  show  the  asymptotic  behavior 

AyR;h“*(jx-0_fc  +  ffo)’.  i  — oo  (4) 

Note  dependence  Xj  =  Ay(r)  and  =  dy(''i<)-  The  independence 
of  the  boundary  conditions  Horn  r  meant  d/,r  satisAes  (3). 


((9?  +  9l».<S,)  =  -Vr- 

DeAne  the  following  semi-inAnite  matrices 


m 


.  A  = 


0  oij 

-oij  0 

-013  “OZ3 


0J3 

0 


(Up 


and  the  diagonal  matrix  A  =  diag(A|,A2,Aj,  )  Uritmg  (10)  and 
(11)  in  matrix  form  gives  the  matrix  equation  for  (I)  h 


Dr-i-ADs  ^AD. 


(13 


The  original  differential  operators  go  to  the  ttansfotmod  mains 
operations  as  -•  +  A  and  +  y  -•  -A.  We  note  the 

following  are  skew-syrrimclric  (Af*  =  is  the  adjoint) 

•b 


(ii“ 


392 


to  the  atagy  cpn*etv»lioa  sUteraenl 
(tf,u)  s  =  3‘B  =  constant  in  r  (15) 


.  t  Th^  calculation  demonstrates  Uiat  mode  coupKos  is  due  et>< 
to  ll>«  »■  dependence  of  ‘Xlic  vertical  operator  5/  +  ?  is 
^iJonned  to  the  diagonal  matrix  A«  vhicb  cannot  contribute 
fj^asse  from  di  to  d>« )  ^  K  via  ofT-dia^onal  elements.  The 
^rTraUse«  hoire%'er,  transf^ms  to  an  r  deri/ative  ploe  an  aaU* 
^Qjnxtric  matrix  A  in  (13).  Now  if  ve  have  pressure  release  top 
pj  a  tip^  bottom,  the  eigenvalues  and  eigenfunctions  are 


-  -  V  o 


mxz 

f^ftrcsio  — 


rde  la^t  b  arbitrar)  )  The  inner  product  in  oj^  in  ($)  is  approxi- 

(l.rfi  M  =  -  >»’)■■: 

p^S  when  /  +  p  +  m  b odd,  and  ey^  s  p  +  j  —  I  when 
j■^p+m  be*eR.  Note  that  the  interaction  couples  a  wide  range 
cf  frequencies,  although  when  various  sums  and  difTerciiees  of  j,  p, 
tad  m  ate  small,  there  is  greater  resonance. 

(D)  If(^>)ist.3.~  decay  fast  enough,  repeated  use  of  (11)  gives 

W;+?1«.W  =  n-V)-  (".#,)  (>6) 

pdjnomial  /’(A)  Comergence  arguments  for  (16)  gi«e  that  a 
tcBctwn  /  may  be  rewritten  in  matrix  form  /(d/+^]u  -•  /(-A)fl, 
vhete 

/(-A)=  d!>s(/(->i).  n-X,).  /(-A,).  ...)  (17) 

A  particularly  impoiUnt  use  of  this  relation  occurs  with  the  fac- 
terlralion  of  (i*d?  +  d*  +  q  +  ^*Iu  =  0  into  the  general  parabolic 
eperators 

(ikdr  -  QKii9r  +  C)  -  rt(d,.<51u  =  0  (1$) 

vhere  Q  =  \/dj +  q+t’  and  (Af,  fi]~  MN  -  NM  is  the  commu- 
Ulor,  Q  transforms  to  the  diagonal  matrix  operator  C  s  —  A 
ttdescribed  by  (17)  Thb  permits  the  exact  calculation  of  the  cor* 
responding  commutator  +  A,CJ  s  C  where 

C^,~  .<*,)  (  )’*. 

G  exists  as  long  as  — A|.  ^  A*.  Thus,  the  ap|>roximation  may  breaV* 
Aown  when  \p  fa  k^,  which  may  be  estimated  by  (4)  In  fact  if  we 
la>e  continuous  spectrum  in  this  range,  there  wilt  necessarily  be  a 
Irealdown  in  the  commutator,  and  possibly  in  the  approximation 
cf  the  eltiplie  equation  by  the  parabolic. 

(C)  Now  use  y/d^  +  A*  -•  VA*  —  A  to  concert 


+  ^\/d?  +  q  +  A*^  u  = 


to  the  matrix  difTerential  equation 

B,,  =  VIi,  Kb  (W) 

i^ck  (he  sign  of  the  radical  in  (20)  to  be  (a)  positive  when  A^  >  Ay 
wd  to  be  (b)  positive  imaginary  when  A*  <  Ay.  Let  Pp  be  the  ma* 
trix  with  ones  on  the  diagonal  in  case  (a),  holds  and  reros  other* 
vise,  and  /«  have  ones  on  the  diagonal  in  case  (b)  and  xeros  other* 
vise.  Because  of  the  choice  of  sign,  Pp  and  P,  are  projections  onto 


the  propagating  modes  and  the  evanescent  modes,  respectively. 

Pp  +  Pt“it  the  identity  matrix.  Write  Bp  s  /^B=:(Ai,...,A«,)^, 
the  propagating  modes  of  B.  For  these  modes  we  may  approximate 
the  radi^ 

(lA)-*  ss  -*  -  (2iA*)-‘A 

fairly  well,  since  ihe  maximum  error  is  j.  Thb  approximation 
gives  (13)  «h«  ve  include  a  complex  exponential  u  =  ue^^'.  In 
case  (b),  however,  (hb  approximation  faib,  since  one  may  show 
(B«  Bt)^  <0  Thus  ut  decays  with  respect  to  v,  as  opposed  to 
energy  conservatioo  in  (15)  Fade  approccimanta  of  the  square  root 
may  also  cause  problems,  in  (hat,  if  ooe  of  the  Ay  is  near  a  pole  of 
the  Fade  approximant,  it  could  induce  rapid  oscillations  rot  due 
to  themginal  operator.  One  way  around  (hb  b  to  use  separate 
approximations  for  the  propagating  and  evanescent  parts.  A  Fade 
app^kxtnutioa  should  work  well  for  the  evanescent  modes[2]. 

(O)  V  in  (20)  may  be  viewed  as  a  generator  of  a  contraction  semi* 
group  of  operatM,  if  it  includes  tome  constant  attenuation.  If  q  b 
a  Marlw  process  of  r,  Ifersh  and  Papanicolaou  {9}  dbcuss  the  exis¬ 
tence  of  an  averaged  operator  (F).  Thb  could  be  used  to  e«timate 
an  average  pirabohc  equation  for  (19). 

4  SUMMARY 

The  lepresentalion  of  the  parabolic  equation  in  range  depen' 
dent  normal  modes  allows  us  to  write  an  infinite  dimensional  ma' 
trix  equation  F^om  (hb  we  see  that  mode  coupling  b  due  to  mode 
interaction  with  the  range  dependence  of  the  index  of  refraction. 
The  generalised  parabolic  equation  (19)  gives  fairly  good  approx* 
imation  of  the  elliptic  wave  equation,  except  in  the  spectral  range 
where  propagating  modes  go  over  to  evanescent  modes  When  one 
approximates  Q  by  simpler  operators,  one  should  approximate  it 
separately  on  the  propagating  and  evanescent  parts  of  the  spec¬ 
trum,  and  (hen  mtroduce  the  coupling  between  these  parts. 
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